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Hypobaric  Decompression  Sickness  Model  Development  (Part  I); 
Diffusion  of  Inert  Gas  from  a  Viscoelastic  Fluid  (Blood) 
into  an  Expanding  Gas  Phase 


Introduction 


In  1989,  an  altitude  decompression  model  development  program  was  initiated  at  the  Armstrong 
Laboratory.  The  immediate  goal  of  this  effort  is  to  define  the  architecture  (or  framework)  and  the  algorithms 
(or  software)  for  the  decompression  model.  The  final  product  will  be  a  software  package  that  will  serve  as 
the  altitude  decompression  "standard*  and  will  lead  to  hardware  development  to  provide  the  USAF  with 
DCS  risk  assessment  capability  for  a  variety  of  operational  settings.  This  report  is  one  in  a  series  of  several 
steps  taken  towards  the  completion  of  this  goal. 

In  ct'jdies  on  cavitation  and  cavitation  damage,  it  is  necessary  to  use  mathematical  expressions 
to  characterize  pressure  and  velocity  fields  in  the  neighborhood  of  a  growing  or  collapsing  gas-filled  cavity 
in  a  fluid.  Our  emphasis  is  on  the  problem  of  spherical  bubbles  in  an  incompressible  viscoelastic  liquid. 
Attempts  have  been  made  recently  to  associate  the  number  of  bubbles  and  their  size,  with  the  existence 
and  severity  of  decompression  sickness  (DCS).  It  is  a  fact  that  bubbles  appear  when  dissolved  gases  are 
supersaturated  in  the  blood  and  the  tissue  during  DCS.  Bubbles  have  b^n  seen  with  various  methods 
(Ultrasound  imaging  techniques),  and  It  is  commonly  believed  that  their  size  and  numbers  are  involved  in 
DCS.  The  direct  relationship  is  yet  to  be  discovert  because  it  is  still  not  known  how  to  relate  critical 
variables  (that  will  be  able  to  pre-inform  us  if  and  when  DCS  will  occur  and  where)  with  the  actual 
physiological  processes.  The  details  of  the  processes  underlying  some  of  the  manifestations  of  altitude 
decompression  sickness  remain  unknown,  but  the  basic  mechanism  is  undoubtedly  the  supersaturation  of 
the  tissues  and  blood  with  nitrogen.  DCS  may  arise  either  in  flight  or  during  exposure  to  reduced 
atmospheric  pressure  in  a  hypobaric  chamber.  The  dynamic  representation  of  the  diffusion  process,  and 
the  understanding  of  its  solutions,  analytical  (if  analytical  solutions  are  possible)  or  numerical,  will  provide 
an  important  tool  to  better  understand  this  complex  physiological  system. 


Effects  of  viscoelasticity  and  compressibility  on  the  bubble  dynamics  have  been  previously 
investigated.  Barlow  and  Langlois'^’  considered  the  diffusion-fed  growth  of  a  spherical  gas  buU>le  into  a 
Newtonian  viscous  liquid  under  isothermal  conditions.  Street'^  performed  a  theoretical  study  of  bubble 
growth  in  a  viscoelastic  medium  (represented  by  the  Oldroyd  three-constant  model)  by  considering  the 
momentum  transfer  process,  but  neglected  a  diffusion  process  between  the  gas  and  the  liquid  phase. 
Tanasawa  and  Yang'^  also  investigated  the  collapse  of  a  gas  bubble  suspended  in  a  viscoelastic  liquid 
represented  by  the  Oldroyd  three-constant  model.  Like  Street,  they  neglected  the  diffusion  process  between 
the  gas  and  the  liquid  phase  because  they  made  the  assumption  that  the  gas  inside  the  bubble  was 
undergoing  a  reversible  polytropic  process  (the  instantaneous  gas  pressure  inside  the  bubble  was  a 
function  of  a  polytropic  exponent  and  not  coupled  with  a  diffusion  process).  By  taking  into  account  the 


thermodynamic  behavior  of  the  gas  phase  inside  the  bubble,  they  observed  that,  when  the  bubble  _ 
oscillation  occurs  during  collapse,  viscous  damping  is  less  important  in  viscoelastic  liquids  than  it  is  in 
purely  viscous  liquids.  They  obsenred  further  that  a  gas  bubble  collapses  faster  under  adiabatic  conditions  ' 
than  under  isothermal  conditions.  We  will  investigate  a  similar  case,  but  with  the  exception  that  we  will  take 
under  consideration  the  fact  that  the  blood  diffusion  rate  cannot  be  characterized  either  fast  nor  stow  (for 
this  reason  the  diffusion  equation  cannot  be  neglected).  If  the  diffusion  coefficient  of  a  gas  dissolved  in  the 
suspending  medium  lies  in  the  intermediate  range,  we  must  solve  the  hydrodynamic  and  the  diffusion :: 
equation  simultaneously. 
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Assumptions 


We  assume,  first  of  all,  that  a  single  bubble  of  gas  is  growing  in  an  otherwise  unlimited  volume  of 
liquid  in  the  absence  of  foam  dynamics  (in  other  words  we  assume  that  bubbles  won’t  compete  for  space 
with  others  growing  in  close  proximity).  We  assume  that  the  liquid  is  a  viscoelastic  fluid,  and  we  will  use 
an  Oldroyd  type  B  model  to  represent  the  rheological  properties  of  the  blood.  We  also  assume  that  the 
liquid  is  incompressible  (low  velocity  gas  flow  can  be  assumed  to  be  incompressible);  the  gas  is  ideal,  it 
diffuses  according  to  Pick’s  law.  The  concentratbn  inside  the  bubble  is  related  to  the  concentration  just 
outside  the  bubble  according  to  Henry’s  law,  and  only  a  negligible  amount  of  gas  is  absorbed  on  the  bubble 
wall.  Once  we  reach  the  desirable  altitude,  the  initial  radius  value  of  the  gas  pocket  will  change  according 
to  Boyle’s  law  (This  change  occurs  once  and  is  discrete).  The  Laplace  relationship  is  valid  because  we 
emphasize  growth  and  not  nucleation.  The  bubbles  are  assumed  to  be  spherical. 

Hydrodynamics  of  the  Model 

First  we  must  define  the  hydrodynamics  of  the  problem.  When  a  bubble  grows  in  a  liquid,  there  will 
be  a  velocity  field  within  the  liquid  at  the  gas-liquid  interface,  which  in  turn,  generates  a  stress  field  tending 
to  retard  the  bubble’s  growth.  Basically,  when  the  fluid  is  being  pushed  by  the  growing  bubble,  it  creates 
a  velocity  field,  which  is  equal  to  the  rate  of  change  of  the  bubble  radius.  The  spherical  symmetry  of  the 
situation  (the  bubbles  are  spherical),  makes  H  convenient  to  choose  a  spherical  coordinate  system  with  its 
origin  at  the  bubble  center.  The  relationship  of  our  hydrodynamic  equation  and  the  diffusion  equation  is 
governed  by  the  diffusion  rate.  Peticolas*^  calculated  the  limiting  cases  of  extremely  rapid  and  extremely 
slow  gas  diffusion,  and  he  came  to  the  following  conclusion;  if  the  diffusion  is  sufficiently  rapid,  the  bubble 
pressure  remains  constant  and  the  growth  is  determined  by  the  hydrodynamic  equations  alone.  Here  the 
treatment  of  the  moving  boundary  is  important.  On  the  other  hand,  if  the  diffusion  is  slow,  hydrodynamic 
effects  become  negligible,  and  the  diffusion  equation  alone  describes  the  process  fairly  accurately.  An 
asymptotic  solution  of  the  second  case,  for  viscous  liquids,  was  derived  by  Barlow  and  Langlois'^’  in  1962. 
We  will  solve  the  intermediate  case  where  the  hydrodynamic  equation  is  coupled  with  the  diffusion 
equation. 

The  Model 


The  velocity  field  generated  in  the  liquid  can  be  described  completely  by  the  Navier-Stokes 
equations.  If  we  assume  that  it  has  only  a  radial  corrponent  v(r,t),  it  reduces  to:  v  =v,(r,t),  Vg=v,=0. 

From  the  conservation  of  momentum  equation: 


P 

V 


dt  'dr 


dr  ^  ar' 


(1) 


From  the  continuity  equation  in  spherical  coordinates  (See  appendix  1); 


(2) 


Integrating  equation  (2)  from  R(t)  to  r  we  get: 


(3) 
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Notice  that  we  treat  a  stationary  bubble,  with  respect  to  the  blood  flow  (creeping  flow  is  considered  in  the 
veins).  We  can  assume,  because  of  the  biomolecular  structure  of  the  fluid,  that  the  bubble  (like  a 
suspension)  is  just  hitching  a  ride  with  the  blood.  Once  it  is  nucleated,  it  doesn’t  move  independently.  It  just 
follows  the  flow  (Buoyancy  forces,  gravity  effect  are  omitted),  and  it  grows  in  an  infinite  mass  of 
honnogeneous  incompressible  viscoelastic  liquid.  Incompressibility  and  spherical  symmetry  implies 
T,r+Tee+t,^=0  and  Tgg=x,^.  By  taking  under  consideration  these  two  relationships  and  by  substituting  equation 
2  and  3  in  equation  1 ,  we  integrate  equation  1  (with  respect  to  r)  from  R  to  infinity.  The  result  is  the 
following  integrodifferential  equation; 


r 

=PX^-P(“) 


(4) 


The  above  equation  represents  the  pressure  distributbn  in  the  blood.  The  dynamic  equation  governing  the 
growth  or  collapse  of  the  bubble  may  be  obtained  by  replacing  r  with  R.  The  balance  of  forces  at  the 
bubble  surface  (r=R)  requires  that; 


(5) 


in  which  the  radial  normal  stress  acting  on  the  bubble  surface  due  to  the  gas  phase  viscocity  is  neglected. 
Similarly,  the  balance  of  forces  at  r=oo  will  be;  -p(‘«)+T„(«>)=-p„  (Relationship  1 ). 

In  the  present  study,  the  three-parameter,  Oldroyd  model  is  employed  to  represent  the  rheological 
behavior  of  the  blood.**’  The  B  fluid  model  is  proposed  here  as  the  simplest  case  that  might  apply  to  blood 
(shear  thinning,  positive  Weissenberg  effect).  This  model  allows  a  reduction  to  a  three-constant  model  but 
results  in  both  the  elimination  of  the  extra  normal  stress  terms  and  the  modification  of  the  term  ige  in  a 
viscometric  flow.  However,  since  it  does  not  appear  that  one  can  determine  the  additional  constants 
experimentally,  the  model  is  employed  as  a  reasonable  approach  in  the  spirit  of  most  nrxxiels  proposed  for 
the  complex  blood  fluid.  Other  more  recent  models  for  the  blood  are  the  three-dimensional  dyadic  Walburn- 
Schneck  constitutive  equation  derived  by  Easthope®  and  the  non-linear  Maxwell  model  by  Quemada'^^. 
Both  of  these  models  will  be  examined  in  the  future. 


The  full  constitutive  equation  for  the  blood  was  given  by  Phillips  and  Deutsch  in  their  1 975  paper.'*’ 
Here  we  will  treat  the  three-parameter  model,  which  can  be  written  in  spherical  coordinates; 


(6) 


where  D/Dt  is  the  material  derivative  defined; 


<7, 

Ti„  means  that  the  stress  variable  is  a  tensor,  and  it  is  in  a  matrix  form.  Because  of  our  assumptions,  x  is 
equal  to  the  following  matrix; 
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(8) 
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is  the  characteristic  stress-relaxation  time,  is  the  shear  viscosity,  is  the  characteristic  strain- 
relaxation  time  (X,>X2>0),  and  d|^  is  the  rate  of  strain  tensor  which  can  be  expressed; 


drr  0  O' 
0^0 
Q  Q  d 


w 

dr 

0 

0 


0  0 


r) 


(9) 


The  special  case  in  which  X,=X2=0  corresponds  to  a  Newrtonian  fluid.  Remember  equation  6  is  in  matrix 
form.  The  stress  tensor  is  being  represented  by  equatbn  8;  consequently,  equation  6  is  actually  a  system 
of  3  differential  equations.  Because  we  are  talking  about  a  velocity  field  in  the  liquid,  we  are  interested  only 
in  the  radial  stress  component.  Thus,  the  matrix  form  of  equation  6  can  be  reduced  to  the  following: 


(10) 


We  use  the  same  technique  that  Tanasawa  and  Yang  used  in  their  1970  paper^^,  in  order  to  avoid 
nonlinear  convective  terms  by  introducing  a  new  independent  variable:  y=(r’-[R(t)f)/3.  The  material 
derivative  defined  in  equation  7  in  the  (r,t)  coordinate  system  may  be  reduced  to  the  following  form  in  the 
(y,t)  coordinate  system: 


D_d  ^  f  ^  ^ 
dy  dydt 


Therefore  equation  10  can  be  written: 


dt 


+'c^,0=-4ti^Z+A.24 


(11) 


(12) 


Where  2=dZ/dt  and 


(Sy.PP) 


(13) 
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Thus,  through  this  transformation,  the  non-linear  convective  terms  in  equation  10  are  eliminated.  The 
normal  stress  component  can  be  obtained  by  integrating  equation  12  as; 


3y*PP 


di. 


(14) 


Finally  equations  4,5,14  and  relationship  1  are  combined  and  the  resulting  equation  is  integrated  with 
respect  to  y.  It  yields; 


(15) 


i2tio  f  ^ 


This  relationship  is  an  integro-drfferential  equation  that  tells  us  how  the  bubble  radius  varies  with  respect 
to  time.  In  addition  to  R,  we  have  another  variable  whose  rate  we  also  have  to  define.  The  pressure  p,  at 
any  point  in  the  liquid  is  a  function  of  r  and  t.  However  the  pressure  p^  (pressure  of  the  ideal  gas  forming 
the  bubble)  is  only  a  function  of  time.  We  need  an  expression  that  relates  Pg  to  bubble  radius  R.  Such  a 
relationship  can  be  derived  from  a  mass  balance  performed  on  the  gas  dissolved  in  the  liquid  phase  by 
Pick’s  law  of  diffusion.  Before  we  proceed  in  describing  the  process,  it  should  be  mentioned  that  the  bubble 
pressure  Pg(t)  is  related  to  the  amount  of  gas  within  the  bubble  at  time  t  according  to  the  ideal  gas  law.  If 
we  assume  isothermal  conditions,  (in  our  case  this  assumption  is  perfectly  valid,  because  we  keep  the 
cociqait  at  constant  temperature);  Pg(t)=Apg(t)(Relationship  2),  where  p  is  the  concentration  (density)  of  the 
gas  forming  the  bubble,  and  where  A  is  a  constant  given  by  A=RgT^.  Constant  Rg  is  the  universal  gas 
constant,  T  is  the  absolute  temperature,  and  M  is  the  nwlecular  weight  of  the  gas  (nitrogen  in  our  case). 

Assume  that  at  time  t=0,  a  homogeneous  concentration  of  gas  is  dissolved  throughout  the  liquid. 
As  mentioned  above,  the  concentration  C  obeys  Pick’s  law  of  diffusion,  and  in  spherical  coordinates  has 
the  form; 


where  V=R^(dR/dt)/r^  and  D  is  the  diffusion  constant.  The  material  derivative  is  used  to  account  for 
convection  of  gas  by  the  moving  liquid.  At  the  time  the  pitot  first  reaches  altitude  the  bubble  (actually  a 
nuclei,  previously  in  equilibrium  at  sea  level  conditions,  being  created  through  heterogeneous  nucleation 
or  tribonucleation)  will  grow  instantly  (Boyle’s  law),  resulting  in  a  discrete  drop  in  pressure  inside  the 
bubble.  As  the  gas  bubble  grows,  the  pressure  (concentration)  inside  the  bubble  will  decrease,  giving  rise 
to  a  ’concentration  gradient  at  the  gas-liquid  interface,  which  in  turn,  affects  the  mass  transfer  process 
(This  concentration  gradient  depends  on  the  altitude  that  the  pitot  was  exposed  to,  and  the  prebreathe^ 


’t  prebreathe=100%  Oxygen  intake  prior  to  exposure  to  maximum  altitude. 
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protocol).  Hence  equation  16  must  be  solved  simultaneously  with  an  expression  that  relates  the  solute 
concentration  just  outside  the  bubble  wall,  C,.(t),  to  the  concentration  inside  the  gas  bubble,  Cg(t).  Henry’s 
law  will  provide  that  relationship  for  us:  C|(R,t)-C  Jt)=KpPg(t)=KpApg(t).  Then  pg(t)=C^(t).0<pA(Relationship  3). 
Constant  K,.  is  the  Henry’s  law  constant.  The  initial  condrtions  tor  R{t)  and  C,(r,t)  may  be  written  as  R(0)=R„, 
dR(0)/dt=0,  and  C,(r,0)=C„. 

When  the  arterial  nitrogen  concentration  drops,  either  because  of  altitude  exposure  or  because  of 
an  increase  of  the  oxygen  concentration  in  the  breathing  air,  the  venule  nitrogen  partial  pressure  depends 
on  the  capillary  diffusion-perfusion  gas  exchange  process  between  the  blood  and  the  tissue.  This  exchange 
process  serves  as  a  control  mechanism  of  the  venule  nitrogen  partial  pressure.  An  increase  in  the  venule 
nitrogen  partial  pressure  results  in  a  concentration  gradient  between  the  blood  and  the  bubble,  and  the 
activation  of  the  growth  process.  For  this  reason,  our  analysis  is  better  applied  in  bubbles  that  nucleate  and 
grow  in  the  veins.  Here,  the  nitrogen  concentration  does  not  vary  greatly  with  distance;  that’s  why  we  will 
assume  that,  at  a  relative  distance  from  the  bubble,  the  concentration  of  nitrogen  will  be  constant.  (In 
general  this  is  not  true.  In  order  to  understand  how  the  nitrogen  concentration  changes  axially,  radially  and 
due  to  convection  in  the  veins,  we  must  fully  solve  and  understand  the  parabolic  transport  model, 
describing  capillary  gas  exchange  properties.  This  has  been  the  topic  of  several  papers  by  Goresky*'®’,  and 
Bassingthwaighte*’'**).  Translating  this  assumption  into  mathematics  provides  us  with  our  first  boundary 
condition; 


X\mCI/.ti=C„  (17) 


At  a  sufficiently  great  distance  from  the  gas-liquid  interface,  the  effect  of  the  growing  bubble  should  be 
negligible;  thus,  the  concentration  at  any  time  is  nearly  equal  to  the  initial  concentration.  The  second 
boundary  condition  representing  the  mass  flux  at  the  boundary  can  be  expressed; 


In  order  to  gain  information  on  bubble  growth  rate,  we  must  solve  equations  15,  16  and  18 
simultaneously.  In  so  doing,  we  shall  adopt  the  integral  or  moment  method  (one  member  of  a  class  of  so- 
called  "weighted  residual"  methods)  introduced  by  Rosner  and  Epstein'^',  which  permits  us  to  obtain  an 
expression  for  the  bubble  radius-time  relationship  in  closed  form.  It  provides  a  practical  and  sufficiently 
accurate  approach  to  these  nonlinear,  multiparameter  boundary  value  problems.  As  Rosner  and  Epstein 
mentioned  in  their  paper,  this  method  is  p  dicularly  well  suited  to  problems  in  which  the  dependent  variable 
(solute  concentration  in  our  case)  at  the  liquid-vapor  interface  is  time-dependent.  Applying  this  method,  we 
manage  to  reduce  our  partial  differential  system  to  a  system  of  ordinary  differential  equations  in  which  time 
is  the  (only)  independent  variable. 

Multiplying  both  sides  of  equation  16  by  and  integrating  the  resulting  equation  with  respect  to  r 
from  r=R  to  r=R+5  yields: 


(19) 
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in  which  8  is  a  thin  concentration  boundary  layer.  Next  we  integrate  each  term  taking  under  consideration 
the  foibwing  relationships:  r^,=R^R,  C|(R+8)=Co,  C|(R)=C^,  dC/dr=0  at  r=R+8.  Thus: 


(20) 


where  is  the  concentratbn  of  the  gas  dissolved  in  the  blood  far  from  the  bubble,  whbh  as  we  mentbned 
previously,  may  be  considered  to  remain  constant  during  the  entire  period  of  bubble  growth.  Combining 
equations  18  and  20,  and  integrating  the  resulting  equation  with  respect  to  time: 


c,-c. 


(21) 


where  Pg„  is  the  density  of  the  gas  bubble  at  t=0,  arb  R„  is  the  initial  bubble  radius.  Now  consider  the 
foibwing  class  of  profiles  for  the  normalized  solute  concentration  defect: 


_  (1 

C  -C  Q~for-r>R*b 


(22) 


in  which  C^=C^  because  of  boundary  conditions,  and  C=(r-R)/8.  Substitutbn  of  equation  22  into  the  left- 
hand-side  of  21  yiebs: 


Cc-C, 


(23) 


We  assume  that  #VRc1 ;  thus,  we  can  disregard  the  terms  (8/R)^  and  (8/R)®  in  equatbn  23  because  they 
are  very  small.  Then,  equatbn  23  becomes: 


Substituting  equation  22  into  the  right  hand  of  equatbn  18  and  eliminating  5  from  the  resulting  expressbn, 
with  the  aid  of  equatbn  24: 


6p?D{C„-Cji^R* 


(25) 


Equatbns  15  and  25  must  be  solved  simultaneously  for  R(t)  and  C^(t)  with  the  aid  of  relationships  2  and 
3,  using  the  foibwing  initial  conditbns:  R(0)=R„,  dR(0)/dt=0,  CJ0)=C„. 
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Discussion 


We  have  developed  a  bubble  growth  model  for  the  blood.  Special  attention  must  be  paid  to  two 
v'*  ihe  inrtial  values  involved  in  our  computations.  The  first  one  is  the  dR(0)/dt  value.  We  assumed  that  this 
value  was  equal  to  zero,  but  this  assumption  was  not  totally  valid.  If  we  look  into  the  nucleation  process 
that  takes  place  before  growth,  the  size  of  the  embryo  (nucleus)  constantly  changes  until  it  reaches  a 
critical  radius  size  Embryos  that  never  reach  this  size  soon  collapse.  The  ones  that  do  manage  to  reach 
this  r-critical  radius  size  (through  energy  fluctuations),  are  destined  to  grow.  In  our  analysis,  we  “peak-up" 
the  growth  process  exactly  when  the  embryo  crosses  the  r-critical  borderline  and  becomes  a  bubble.  This 
threshold  means  that  the  rate  of  radius  gro\wth  at  t=0  (which  is  the  starting  time  for  the  growth  process,  but 
not  for  the  nucleation  process)  is  not  actually  zero.  However,  there  is  a  case  in  which  the  initial  radius  rate 
of  change  might,  in  fact,  be  zero.  There  are  theories  that  claim  the  pre-existence  of  nuclei  in  stable  form 
In  either  case,  the  variable  initial  radius  growth  rates  do  not  affect  our  final  radius  computations.  Barlow 
et  al.’’*  solved  analytically  the  initial  stage  of  growth,  at  a  low  Raynolds  number,  for  the  viscous  case,  and 
he  concluded  the  following:"  Except  for  the  first  moments  of  growth,  the  solution  is  virtually  independent 
of  the  value  chosen  for  dR(0)/dt".  Thus,  none  of  the  two  above  assumptions  will  produce  a  significant  error 
in  our  calculations.  Nevertheless,  in  this  study,  we  will  assume  that  dR(0)/dt  is  equal  to  zero. 

The  second  initial  value  that  we  have  to  pay  extra  attention  to  is  the  initial  nitrogen  concentration 
in  the  blood  (on  the  venous  side).  This  issue  is  slightly  more  complicated  than  the  previous  one. 
Prebreathing  oxygen  before  or  during  the  flight  (and/or  exposure  to  altitude)  affects  the  pilot’s  blood  and 
tissue  nitrogen  concentration  levels.  All  body  tissues  in  equilibrium  wKh  ambient  nitrogen  will  respond  to 
increased  or  decreased  nitrogen  partial  pressure  in  the  breathing  medium  by  absorbing  or  eliminating 
dissolved  gas  in  the  tissues.  The  nitrogen  partial  pressure  of  the  blood  in  the  arterial  side  can  decrease 
very  rapidly;  in  the  case  of  100%  Oj  prebreathe,  it  can  drop  to  zero,  creating  a  large  concentration  gradient 
between  the  tissue  and  the  blood  in  the  capillary  level.  Consequently,  the  blood  begins  to  wash  out  nitrogen 
from  the  tissue  through  diffusion.  The  tissue  can  be  characterized  by  a  spectrum  of  theoretical  tissue  types 
that  represent  "fast"  and  "stow"  tissues.  The  outgoing  nitrogen  partial  pressure  in  the  blood  (outgoing 
meaning  from  the  capillary  level  to  the  venous  side,  basically  the  venous  nitrogen  partial  pressure),  can 
vary  depending  on  the  "half-time"  of  the  tissue.  A  theoretical  tissue  type  is  characterized  by  the  tissue  half¬ 
time;  the  time  required  for  a  tissue  to  respond  to  a  change  in  ambient  nitrogen  partial  pressure  by  giving 
off  or  absorbing  nitrogen  until  the  initial  difference  between  the  tissue  nitrogen  partial  pressure  and  ambient 
nitrogen  partial  pressure  is  reduced  by  one-half.  Half-time  tissues  are  very  helpfu!  if  one  follows  the 
exponential  function  approach.  (Since  we  haven’t  developed  our  own  parabolic  transport  model  yet,  we  will 
use  the  exponential  function  approach  to  calculate  the  nitrogen  concentration  near  the  venule  walls.) 

The  outgoing  nitrogen  partial  pressure  must  be  quantified  before  we  can  resume  our  computations. 
The  nitrogen  gas  dynamics  and  exchange  in  the  capillaries  is  not  a  simple  mathematical  interpretation.  For 
the  purpose  of  this  report,  we  will  compute  bubble  growth  at  the  venule  walls  because  that’s  where  the 
nitrogen  concentration  is  at  its  maximum  (maximum  concentration  in  the  blood  results  in  a  maximum 
concentration  gradient  between  the  blood  and  the  bubble).  Its  value  will  be  calculated  by  the  use  of  the 
washout  exponential  model  developed  by  Conkin  et  al.'®*  This  method  is  used  because,  even  if  we 
consider  the  perfusion-diffusion  process,  the  boundary  condition  in  the  capillaries  will  be  C(r=a,t)=C,=C„ 
(a=the  radius  value  of  the  capillary;  C,  is  the  tissue  nitrogen  concentration,  near  the  capillary  wall;  C„  is  the 
venule  nitrogen  concentration  near  the  vessel’s  wall  and  close  to  the  venule-capillary  interface).  C,  will  vary 
according  to  the  Conkin  et  al.  equation  and  so  will  the  nitrogen  concentration  in  the  venule  walls  (We  are 
expecting  the  nitrogen  concentration  value  to  be  higher  in  the  venule  walls  no  matter  what  kind  of  analysis 
we  use,  because  of  the  radial  diffusion  process  in  the  capillaries  ).  It  is  logical  (for  the  above  reasott)  to 
expect  most  bubbles  to  form  closer  to  the  walls. 
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Thus,  the  tissue,  depending  on  its  half-time  value,  alters  its  partial  pressure  value  due  to  the 
concentration  gradient  between  the  tissue  and  the  blood.  This  value  can  be  calculated  through  the  standard 
Haldanian  approach  from  the  following  equation; 


PrPo^{P.-P^(^ 


(26) 


where  P,  is  the  nitrogen  partial  pressure  in  the  tissue  after  exposure  for  t  minutes;  is  the  initial  tissue 
nitrogen  partial  pressure;  P3  is  th'?  ambient  nitrogen  partial  pressure  in  the  breathing  medium;  t  is  the 
exposure  time  in  minutes;  and  k  is  the  tissue  nitrogen  partial  pressure  rate  constant  (k=.693/t,,2  where 
is  the  tissue  nitrogen  partial  pressure  half  time).  Now  we  can  solve  our  model  taking  under  consideration 
an  array  of  half-time  tissues. 


Clarifications 

1 .  It  is  customary  to  express  the  gas  concentration  in  a  tissue  in  terms  of  its  partial  pressure  p  and  the 
solubility  a  of  that  gas  in  the  tissue;  C=aP.  In  this  report,  concentration  and  partial  pressure  are  used 
constantly;  hence,  it  was  necessary  to  mention  their  above  relationship  to  avoid  confusion.  Also  we  assume 
that  concentration  changes  in  the  blood  (in  the  veins,  due  to  convection)  are  much  slower  than  bubble 
growth  rates. 

2.  There  was  no  particular  reason  why  we  chose  an  exponential  decay  approach  (a  perfusion  dominated 
approach)  rather  than  Hills’*"'  approach  (a  diffusion  dominated  approach).  Hills’  equivalent  equation  for  pure 
radial  diffusion  car.  be  expressed  as; 


2D,  S„ 

(dPc 

dr. 

(27) 


and  the  above  equation  is  valid  only  if  v  (the  capillary  blood  velocity)  »>  I  (the  capillary  length).  What  each 
of  the  parameters  and  variables  represent  is  not  important  since  it  is  clear  now  that  the  diffusion-perfusion 
model  must  be  solved  in  order  to  thoroughly  understand  capillary  transport  phenomena  and  nitrogen 
concentration  values  in  the  veins.  This  will  be  the  purpose  of  the  second  publication  in  this  series. 

3.  Using  polynomial  profiles  for  the  gas  concentration  is  tantamount  to  assuming  that  the  bubble  grows  in 
an  infinite  pool  of  liquid.  In  situations  where  a  large  number  of  bubbles  nucleate  and  grow  simultaneously 
in  close  proximity,  the  amount  of  liquid  immediately  surrounding  a  bubble  is  finite  and  so  is  the  amount  of 
the  dissolved  gas.  Therefore,  the  hydrodynamics  of  the  growth  and  the  diffusion  process  will  be  different 
from  the  case  of  a  bubble  growing  in  an  infinite  medium.  "The  underlying  assumption  in  adopting  a 
polynomial  profile  to  describe  the  gas  concentration  is  that  at  large  distances  from  the  bubble  interface, 
the  gas  concentration  remains  unchanged  and  is  equal  to  the  initic.!  concentration.  This  assumption  is  not 
valid  for  a  bubble  surrounded  by  a  finite  amount  of  liquid  with  a  limited  concentration  of  dissolved  gas  and 
cannot  predict  the  steady-state  bubble  radius  correctly."*®’  Therefore,  it  is  necessary  to  solve  the  diffusion 
equation  in  its  complete  form  numerically  as  opposed  to  using  similarity  solutions*'®’  or  approximate 
analytical  methods.  In  their  paper,  Michaelides  et  al.*®’  analyzed  the  process  of  mass  diffusion-induced 
growth  of  a  gas  bubble  surrounded  by  a  limited  amount  of  liquid  and  dissolved  gas  and  presented  a 
solution  technique  to  accurately  solve  for  the  concentration  profile  of  the  dissolved  gas.  Bubble  growth 
dynamics  was  predicted  with  this  concentration  profile.  As  mentioned  in  our  assumptions,  we  will  not 
include  foam  dynamics,  so  the  use  of  polynomial  profiles  for  the  gas  concentration  is  sufficient. 
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Model  Summary 


The  equations  that  describe  bubble  growth  in  the  blood,  after  all  the  mathematical  manipulations,  can  be 
summarized  in  the  following  system  of  ordinary  differential  equations: 


12^0 f . 


And  from  the  mass  transfer  process: 


and 

P,(t)=Ap,(t) 

C,(t)=K,P,(t) 

With  the  following  initial  conditions:  R(0)=Ro,  dR(0)/dt=0,  C*(0)=Co 


R  =bubble  radius 
p,  =density  of  the  liquid  phase 
Pg  =density  of  the  gas  inside  the 
bubble 

pgo=density  of  the  gas  bubble  at 
t=0 

Pg  =pressure  inside  the  gas  bubble 
P«,  =pressure  of  the  liquid  phase 
far  away  from  the  bubble  surface 
a  =surface  tension 
A  =ideal  gas  law  constant  (for  a  more 
detailed  definition  look 
at  page  5) 

Kp  =Henry’s  constant 


Notation 

Tip  =shear  viscosity 
A,,  =characteristic  stress 
relaxation  time 
^  =dummy  variable 
t  =time 

Aj  =characteristic  strain 
relaxation  time 
D  =diffusbn  coefficient 
Cp  =initial  concentration 
of  the  gas  dissolved  in 
the  liquid  phase 
C,,  =solute  concentration 
just  outside  the  bubble 
wall 
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Appendix  1 


The  momentum  equation  in  spherical  coordinates  (r.Q,<p)  in  terms  of  x  is  expressed  as  follows  in  the  r 
component: 


dv^ 

'Sr  r  as 

/sine 

dtp  r 

M 

dr 


/sine  de 


(T^ne)+ 


1 

/s/he  dtp 


The  velocity  field  generated  in  the  liquid  will  have  only  a  radial  component  v(r,t),  where  t  is  time  measured 
from  the  instant  of  bubble  formation.  All  the  derivatives  with  respect  to  6  and  (p  will  be  zero.  Also  the 
velocity  components  in  the  0  and  <p  coordinates  will  be  zero.  Taking  under  consideration  the  above 
assumptions,  the  momentum  equation  can  be  rewritten  as  follows; 


a#  '  I  _9  ff' 


dt  'dr)  dr  ^dr 


The  equation  of  continuity  in  spherical  coordinates  ( r.G.cp) 


UCII  I 


As  we  mentioned  before,  all  partial  derivatives  with  respect  to  0  and  (p,  will  be  zero.  The  density  of  the  fluid 
is  not  a  spatial  variable;  therefore,  the  above  equation  can  be  reduced  to  the  following  one: 
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